Scaling theory of driven polymer translocation 
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We present a theoretical argument to derive a scaling law between the mean translocation time 
r and the chain length N for driven polymer translocation. This scaling law explicitly takes into 
account the pore-polymer interactions, which appear as a correction term to asymptotic scaling and 
are responsible for the dominant finite size effects in the process. By eliminating the correction-to- 
scaling term we introduce a rescaled translocation time and show, by employing both the Brownian 
Dynamics Tension Propagation theory [Ikonen et al, Phys. Rev. E 85, 051803 (2012)] and molecular 
dynamics simulations that the rescaled exponent reaches the asymptotic limit in a range of chain 
lengths that is easily accessible to simulations and experiments. The rescaling procedure can also 
be used to quantitatively estimate the magnitude of the pore-polymer interaction from simulations 
or experimental data. Finally, we also consider the case of driven translocation with hydrodynamic 
interactions (His). We show that by augmenting the BDTP theory with His one reaches a good 
agreement between the theory and previous simulation results found in the literature. Our results 
suggest that the scaling relation between r and N is retained even in this case. 



PACS numbers: 

Introduction. The transport of a polymer through 
a nano-sized pore occurs in many biological processes 
and functions, including DNA and RNA translocation 
through nuclear pores, protein transport across mem- 
brane channels and virus injection pQ. The translocation 
process is also envisioned to have several biotechnologi- 
cal applications in gene therapy, drug delivery and rapid 
DNA sequencing [2r[5] . However, despite significant tech- 
nological advances and considerable experimental [21[5J[7] 
and theoretical [SHIP] studies, until very recently the 
connection between fundamental theory and numerical 
simulations and experiments has remained elusive (for 
a review, see, e.g., Ref. j8]). Recently, we proposed a 
theoretical model of driven polymer translocation and 
showed that it explains and reproduces all of the results 
of relevant numerical simulations, hence providing a uni- 
fying view on the physics of driven polymer transloca- 
tion [42j |43] . The model couples a Brownian Dynamics 
(BD) framework for one degree of freedom with the deter- 
ministic Tension Propagation (TP) theory [3J 37 . The 
model has so far been used to study driven transloca- 
tion in the Rouse dynamics picture, where it was shown 
that the pore-polymer interactions induce strong finite 
size effects that persist for extremely long chains up to 
TV « TO 5 g3]. 

The detailed formulation of the BDTP theory is rather 
involved [HI 03] and except for the asymptotic large N 
behavior no explicit analytic solutions to the theory have 
been found. The main purpose of this work is to find 
a scaling type of solution to the BDTP theory. This 
is extremely important in order to obtain an intuitive, 
physical understanding of the driven translocation pro- 



cess. We show here that the pore-polymer interactions 
appear as a correction-to-scaling term to the asymptotic 
behavior of the mean translocation time r as a func- 
tion of the chain length N. We show how this scaling 
relation can be used both to understand and quantify 
the effect of the pore-polymer interactions from molec- 
ular dynamics (MD) simulations or from experimental 
data. Most importantly, by eliminating the correction- 
to-scaling contribution we demonstrate that the corre- 
sponding rescaled exponent reaches its asymptotic limit 
already for chain lengths N < 10 2 - orders of magni- 
tude faster than the conventional exponent. Finally, we 
discuss the effect of hydrodynamic interactions (His) by 
modifying the BDTP model to account for hydrodynam- 
ics in an approximate way. The effect of the His, as calcu- 
lated from the BDTP model, is to shorten the transloca- 
tion time and the corresponding effective exponent due to 
the increased (relative) importance of the pore-polymer 
interactions. The predictions of the BDTP model are 
shown to agree with previous simulation studies found in 
the literature. The result suggests that the scaling solu- 
tion holds even with His, highlighting the importance of 
this result. 

Theory. Although a details of the BDTP theory are 
complex, some central results can be derived from rela- 
tively simple arguments. To find the mean translocation 
time r (or the mean translocation velocity («)), for a 
chain of length TV, one needs to consider the force bal- 
ance between the driving force / and the drag force, as 
suggested, e.g., in Ref. [7]. In principle, the total drag has 
three contributions: the friction due to the cis side sub- 
chain and the solvent, the corresponding friction for the 



trans side subchain and the friction of the chain portion 
inside the pore. For driving forces typically used in exper- 
iments and simulations, the trans side has an almost neg- 
ligible contribution to the overall friction [35J H3] . There- 
fore, the effective friction T can be approximately written 
as r(i) = T} C i S (t) + rj p , where ?7 c i s is the friction of the cis 
side subchain and r\ p is the (effective) pore friction. If the 
length of the pore, l p , is small compared to the contour 
length of the chain, l p <C aN (with a the segment length) , 
the number of segments occupying the pore and thus the 
pore friction r\ p can be regarded as constant during the 
translocation process. 

The friction due to the cis side subchain, on the 
other hand, depends explicitly on the number of mo- 
bile monomers on the cis side. Due to the nonequi- 
librium nature of the problem, this number depends 
not only on the chain length N but also on time 
t [m HH1 IH2 ESI E33 EH HQH43]. Solved as a function 
of time, ^ c i s (t) shows non-monotonic behavior that sug- 
gests the division of the translocation process into two 
stages: the tension propagation stage of increasing fric- 
tion, and the post propagation stage of decreasing fric- 
tion j42l |43|. However, averaged over the whole pro- 
cess, the time-averaged r] cls is approximately given as 
(?7cis) ~ N u , leading to the asymptotic (N — > oo) scaling 
of the mean translocation time asr~ N v+1 |4"0Tl4"5] . 

The full solution of r) C is{t) is rather involved (cf. 
Refs. [HISS]) and does not give much physical insight. 
We will now show that the result (r) ciB ) ~ N v can be ob- 
tained with relatively simple argument, based in part on 
the work of DiMarzio |44j . In equilibrium, the polymer 
assumes a configuration comprising several loops. Before 
the driving force can be transmitted to any given chain 
segment, the preceding loops need to be straighthened. 
Only then can the tension propagate along the chain 
backbone. Therefore, at any given time the friction ?7 c i s 
is due to the motion of the unraveling loop closest to the 
pore, while the rest of the chain is essentially immobile. 
To estimate the average length of the loop, we note that 
for large N, the end-to-end distance of the polymer is 
given by the usual Flory scaling form R ~ aN" . To esti- 
mate the number of times the chain intersects a plane of 
thickness dR parallel to the membrane (see Fig. [I]), we 
note that the number of segments within the plane, dN, 
is proportional to the average line density of monomers, 
N/R, giving dN ~ N x ~ v dR. This is proportional to 
the number of times the chain intersects the plane and 
also approximately the average number of loops (for large 
N). The average length of one loop is thus proportional 
to N/N 1 ^'' — N u . The number of mobile monomers at 
any given time is therefore proportional to N", and the 
friction due to the drag on the cis side is (r/ c is) ~ i]N u , 
where rj is the solvent friction per chain segment. The 
total effective, time-averaged friction is thus of the form 
(r) « Cr\N v + r/p, with C an TV-independent constant. 

For a sufficiently strong driving force /, one can ex- 
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FIG. 1: (Color online) A schematic illustration of the config- 
uration of a polymer tethered to a wall. Intersections of the 
polymer with an imaginary plane of thickness dR oriented 
parallel to the wall are indicated by the red circles. The av- 
erage number of times a polymer of length N intersects the 
plane is proportional to N 1 ^" . 

press the mean translocation time in terms of the av- 
erage translocation velocity (v). Force balance implies 
( v ) = //(r)- Since the driven translocation proceeds 
via gradual uncoiling of the whole chain, the relevant 
length scale of the process is the chain's contour length, 
aN (note that this is in contrast with the assumption of 
Ref. [7j, where uniform contraction of the cis side chain 
suggests the radius of gyration as the relevant length 
scale). Thus, r = aN/(v) = aN(T)/f. With respect 
to the chain length, the complete relation is then 

T^A(f,r 1 )N 1+v + B(f,r 1 )f] p N, (I) 

where A and B are independent of N, and f] p = rjp/rj is 
the dimensionless pore friction. This scaling form is one 
of the main results in this paper. It shows that the influ- 
ence of the pore friction appears as a correction-to-scaling 
to the asymptotic value of the translocation exponent a, 
defined via r ~ N a , which approaches a x = 1 + v from 
below. In addition, Eq. ([I]) suggests that if the dimen- 
sionless pore friction increases (e.g., by either decreasing 
the solvent friction r\ or by reducing the pore radius), 
the exponent a becomes smaller, especially for relatively 
short chains, in agreement with theory and MD sim- 
ulations noi H3 M EU ■ 

Results. To clarify the influence of the pore friction on 
the translocation time, we have solved r time as a func- 
tion of chain length for different pore frictions using the 
BDTP model E3|43]. Using t(N) ~ N a , one can extract 
the effective scaling exponent a(N) as a(N) = 4 j 1 "^- . Be- 
cause of the linear sub-leading term in Eq. Q, a has a 
weak dependence on the chain length, and approaches 
I + v extremely slowly with typical values of the param- 
eters. Naturally, for larger fj p , the dependence on N is 
more pronounced, as shown in Fig. [2] The other parame- 
ter values used in solving the BDTP model were / = 5.0 
(driving force), fceT = 1.2 (temperature), and v — 0.588. 
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FIG. 2: (Color online) Open symbols, main plot: the effec- 
tive exponent a(N) as a function of the chain length N for the 
BDTP model solved for different ratios fj p of pore and solvent 
friction. Solid symbols: the same data rescaled to exclude 
the effect of pore friction, showing collapse of the data to the 
asymptotic exponent 1 + v. Inset: the data collapse for differ- 
ent fjp on a (t/N 1+v ,fjpN~ v ) plot used to extract the finite 
size effect of the pore friction (see text for details) . Numerical 
errors are of the order of the symbol sizes, or smaller. 

To quantitatively show that the deviation from the 
asymptotic limit is caused by the pore friction, we sub- 
tract the correction-to-scaling term and define a rescaled 
translocation time as = r — Bfj p N and the corre- 
sponding rescaled exponent as a^(N) = jf^N ■ Accord- 
ing to Eq. ([!]), the rescaled translocation time can be 
expressed as 

rt = (r - Bfj p N) ~ N a \ (2) 

with the exponent cv independent of N at cr « 1 + 
v. In practice, to obtain the rescaled t\ one has to 
find the numerical prefactor B by means of finite size 
scaling. In the inset of Fig. [2] this is done by plotting 
the translocation time data in the form t/N 1+v = A + 
Bfj p N ~ v . The coefficients A and B can then be obtained 
from a simple linear least squares fit. 

The rescaled exponent a* as obtained from the BDTP 
model is shown in Fig. [2] as solid symbols. All the curves 
corresponding to different values of fj p collapse onto a sin- 
gle master curve within the numerical accuracy around 
= 1 + also coinciding with the ideal fj p = solution. 
For very short chains, the collapse is not perfect because 
of secondary finite chain length effects that come into 
play, and also because as a continuum level description 
the BDTP model may not accurately describe very short 
chains. 

We have confirmed that the rescaling scheme works 
for raw data obtained from MD simulations, too, by per- 
forming Langevin thermostatted MD simulations using 
the Kremer-Grest bead-spring model |45) with typical 
parameters found in the literature (for a more detailed 



FIG. 3: (Color online) Open symbols, main plot: the effective 
exponent a(N) as a function of chain length N as obtained 
from MD simulations (for details, see text). Solid symbols: 
the same data rescaled as in Fig. [2j showing collapse to the 
value 1 + v (dashed line). Inset: Finite size scaling plot used 
to extract the contribution of the pore friction. 

description, see, e.g., Refs. [42j|43]). We used the param- 
eters / = 5.0, kgT = 1.2, i] = 0.7, and pore diameter 
dp = 2.0 and length l p — 1.0 (both in units of the seg- 
ment length a), which correspond to the dimensionless 
pore friction of about fj p ps 5.6 [43]. The chain lengths 
in the simulations were 10 < N < 300, with r averaged 
over at least 5000 successful events for each N. 

Plotting the effective exponent a(N) reveals the same 
dependence on N as predicted by the theory. Further- 
more, by performing finite size scaling similar to the 
BDTP case, we see that the rescaled exponent a 1 ' reaches 
the asymptotic value almost immediately; within the sta- 
tistical error, the value is = 1 + v already for N = 20. 
This shows that Eq. ([T]) is valid for remarkably short 
chains and that the dominant finite size effect in the 
driven translocation process is the frictional interaction 
between the pore and the polymer. 

Finally, we would like to point out another important 
feature of the scaling solution. Used in reverse, the fi- 
nite size scaling procedure can be used to estimate fj p if 
its value is unknown for the pore geometry in question. 
By varying fj p so that the rescaled exponent becomes in- 
dependent of N at m 1 + v, the pore friction can be 
measured for any geometry that satisfies l p <C aN . While 
measuring the pore friction directly from the monomer 
waiting time distribution as outlined in Refs. [42, 43 is 
more accurate, the reverse scaling procedure may be used 
even if the waiting time distribution is not available, as 
is typically the case in experiments. 

Because the driven polymer translocation problem is 
inherently a dynamical, non-equilibrium process, it is ex- 
pected that hydrodynamic interactions (His) from sol- 
vent should play a role. In the simplest approximation, 
the His can be included by considering the Zimm type of 
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friction instead of the Rouse friction as in Langevin dy- 
namics simulations. This means writing down the force 
balance condition such that the drag force (and the fric- 
tion ?7 c j s ) is proportional to the linear length of the mobile 
subchain on the cis side, instead of being proportional to 
the number of mobile monomers. The rationale is that 
for intermediate forces, when the cis side subchain adopts 
a shape reminiscent of a trumpet (or a stem-and-flower), 
the solvent inside the trumpet is also set in motion and 
therefore the monomers inside the trumpet do not fully 
contribute to the drag force. Therefore, the overall drag 
force from the cis side subchain should be slightly smaller 
for the Zimm case. However, since the maximum linear 
length of the mobile part of the cis side subchain is still 
given by the end-to-end distance R ~ aN v , we expect 
Eq. (JlJ to hold even in this case, although with a smaller 
prefactor A in the leading term. 

To test this simple argument, we have implemented 
the Zimm friction of the cis side subchain into the BDTP 
model. To this end, one merely modifies the force balance 
equation [e.g., Eq. (Al) in Ref. [12]] so that the the drag 
force is given by the size of the tension blob (£ ) instead 
of the number of monomers within the blob This 
amounts to replacing the force balance equation of Rouse 
type as f£°™ e (x) = i £ x ^{x'Mx^/af^dx' -> 

/drag" 1 ^) = a J-x V v ( x ')dx' , where v(x) is the instan- 
taneous velocity of the monomers at position x, and x = 
—X is the location of the last mobile monomer. One then 
carries out the same numerical implementation as de- 
scribed in Refs. [12 03] using the modified force balance 
equation. The difference between the Zimm and Rouse 
frictions in the context of the tension propagation formal- 
ism is also further discussed in Refs. [3U |3j>J 133 EH HI] • 
The comparison between the Rouse and Zimm dynam- 
ics is shown in Fig. [4] In both cases the BDTP model was 
solved with / = 5.0, k B T = 1.2, 77 = 0.7, and fj p « 7.14 
for up to N = 5 • 10 9 . For the Rouse friction, the effective 
exponent a(N) is consistently larger than for the Zimm 
friction, and eventually approaches 1 + v for large N. 
For the Zimm friction, the approach to the symptomatic 
limit is considerably slower. In fact, even for N ss 10 9 , 
the numerical value (a s» 1.57) still increases with N. To 
illustrate the importance of the pore friction term to the 
slow convergence, we have solved the model with zero 
pore friction (fj p — 0), showing considerably faster con- 
vergence. Altogether, the numerical results in the large 
N limit seem to be in agreement with 1 + v, which has 
also been predicted to be the asymptotic value using an- 
alytical approximations to the tension propagation the- 
ory ma hx]. 

In the short chain regime, we obtain a(N — 100) f=s 
1.37 ± 0.01 and a(N = 100) « 1.31 ± 0.01 for the 
Rouse and Zimm cases, respectively. Although a de- 
tailed comparison with experiments or hydrodynamical 
simulations is difficult due to the lack of knowledge 
on the pore friction (fj p ), the numerical value of 1.31 
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FIG. 4: (Color online) Main plot: The effective exponent 
a(N) as a function of chain length N for the BDTP model 
with Rouse friction (black squares) and Zimm friction (blue 
circles) for rj p ~ 7.14, and for Zimm friction and fj p = (light 
blue solid line). Inset: The mean translocation time r as 
a function of chain length N for the Rouse (solid black line) 
and Zimm friction (blue dashed line) for rj p ~ 7.14. Numerical 
errors are of the order of the symbol sizes, or smaller. 



seems to be in good agreement with the experimental 
(a « 1.27 ± 0.03, Ref. [7]) and lattice-Boltzmann simu- 
lation results (a ps 1.28 ± 0.01, Ref. [35]). In particular, 
the difference in the exponents measured with and with- 
out His (0.06 ± 0.02) matches the difference reported in 
Ref. [33] (0.08 ± 0.04). In addition, in agreement with 
computer simulations [18, 33], the overall translocation 
time is reduced by the addition of the His, as shown in 
the inset of Fig. [4] 

Conclusions. In this work, we have proposed using 
theoretical scaling arguments that the mean transloca- 
tion time of a polymer driven through a nanopore can be 
written in the form r w A(f, rj)N 1+ " + B(f, T])fj p N. The 
first term is derived from the out-of-equilibrium dynamics 
of the cis side subchain and dominates for large N, while 
the second term stems from the interactions between the 
polymer and the pore and remains significant for the 
typical chain lengths in both experiments and computer 
simulations. This unified scaling form is an important 
physical result and powerful tool for analysis of driven 
translocation. By eliminating the correction-to-scaling 
term one can isolate and quantify the effect of pore fric- 
tion by means of finite size scaling. We have demon- 
strated by using both a theoretical model of translocation 
dynamics and molecular dynamics simulations that the 
rescaled exponent reaches the asymptotic limit already 
for extremely short chains (N < 100), whereas the con- 
ventionally defined exponent does not. In addition, we 
argue that in the presence of hydrodynamic interactions, 
the translocation time becomes shorter but can still be 
expressed as a sum of the two terms. We present results 
from the theoretical model proposed in Refs. 42, 43] with 
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hydrodynamic interactions, obtaining quantitative agree- 
ment in the scaling exponent a with both theoretical and 
experimental results reported in the literature. 
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